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Abstract. The permanent of a square matrix is defined in a way similar to the determinant, but without 
using signs. The exact computation of the permanent is hard, but there are Monte-Carlo algorithms that 
can estimate general permanents. Given a planar diagram of a link L with n crossings, we define a 7?! X 7n 
matrix whose permanent equals to the Jones polynomial of L. This result accompanied with recent work of 
Freedman, Kitaev, Larson and Wang ( IFKLWl ) provides a Monte-Carlo algorithm to any decision problem 
belonging to the class BQP, i.e. such that it can be computed with bounded error in polynomial time using 
quantum resources. 



1. Introduction 

The permanent of an n x n matrix A — (aij ) is defined to be 

n 

per(A) = ^ J|a,^(,) 

where Sym„ is the permutation group on {f , . . . , n}. It is well known that if A is the (V^ x W), 0, 1 adjacency 
matrix of a bipartite graph G = {V,W,E), then per(A) is the number of perfect matchings of G. The 
permanent of A is syntactically similar to its determinant det{A), which is a signed variation of the above 
sum. This mild sign variation leads to a radical change in computability: computing permanents is hard 
fsee ll.3p . whereas determinants can be computed in a polynomial time. 

Here we are interested in the Jones polynomial ([Jlj) which is a celebrated invariant of links in 5'^. A 
link is a disjoint union of embedded circles in 3-space. The Jones polynomial J of a link can be uniquely 
characterized by the following skein relation: 

„2 7 / ^ / \ „-2 7 / \ / \ _ /„ „-l^ 



together with the initial condition J(unknot) = q + q^^ ■ It follows that the Jones polynomial of a link can 
be computed in exponential time (with respect to the number of crossings). 

The Jones polynomial has a statistical mechanics definition, as a state sum which uses as input a planar 
diagram of a link. Our aim is to give a permanent formula of the Jones polynomial. 

Consider a diagram of an oriented link L, that is an oriented, 4-valent plane graph, where each vertex 
has a crossing structure of one of two types: 



positive: / ; or negative: 





We form a graph Dl , from the link diagram , by replacing a neighbourhood of each crossing of with 
one of the graphs shown in Figure [TJ We may say that Dl is a blown- up version of D^. is an immersed 
directed graph with 7n vertices, where n is the number of crossings of D^. We will refer to the the two 
graphs shown in Figure [H minus the incoming and outgoing edges (which come from the link diagram), as 
gadgets. We emphasize the fact that the four edges in Dl that enter and exit a gadget are all parallel when 
they meet the gadget. 
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In Definition 12.21 below, we define local weights on the edges of Dl- Let Ml denote the adjacency matrix 
of weights oi Dl. has size the number of vertices of Z?^, and the entry of is the weight of the 
corresponding directed edge (ij) oi Dl- 




Figure 1. Gadgets at positive and negative crossings. 

Next theorem contains the main result of the paper. An inspiration comes from the observation (see jGL| . 
|BG| ) that the weight system associated with the colored Jones function is a permanent. 

Theorem 1. For every link L we have 

J{L) = q-^^^^^^q'°'<^^^pei{ML). 

1.1. Acknowledgement. The first author wants to thank D. Aharonov, A. Barvinok, S. Basu, P. Bucr- 
gisser, S. Garoufalidis, M. Jcrrum, M. Kiwi and P. Tetali for enlightening discussions. 

1.2. Jones and quantum computing. We follow an exposition of the results of Freedman, Kitaev, Larsen 
and Wang (see |FKLW| ) written by Bordewich, Freedman, Lovasz and Welsh (see in particular proof of 
Theorem 5.1 in |BFLW| ). 

Suppose that we have a BQP language and an input x. Then we can construct a link L, of size polynomial 
in \x\, such that if x is in the language then | J(L, e^'^*/^)| < [2]g^'~''^0.39. On the other hand, if x is not in 
the language then | J(L, e^'^*/^)| > [2]g'^'^"^0.65. We use the standard notation [2]5 — 2cos(7r/5). 

It is also shown in [FKLW| that this kind of approximation of the Jones polynomial (called the additive 
approximation in (BFLWj l is in BQP. This was proven in a different way also by Aharonov, Jones and Landau 
[A JLj . The BQP-hardness of approximating the Jones polynomial is also discussed by Aharonov and Arad 
in [Ml- 
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1.3. Exact computation of Jones polynomial and permanent. The complexity class #P consists of 
the counting versions of the decision problems in NP; an example of a problem in #P is: given a graph, 
how many Hamiltonian cycles it has? It is considered very unlikely that = P. Exact computation of the 
permanent of 0, 1 matrices is #P-complete (see Valiant |V1) ) and exact computation of the Jones polynomial 
J(L, t) is #P— hard except when t is a root of unity of order r € {1, 2, 3, 4, 6} (see Jaeger, Vertigan and Welsh 
[ J VW| ) . Kuperberg (see |Kup| ) showed that stronger than additive approximations of the Jones polynomial 
are hard. This contrasts with the approximations of the permanents, as we explain next. 

1.4. Monte-Carlo algorithms for a permanent. Jerrum, Sinclair and Vigoda constructed fully polyno- 
mial randomized approximation scheme (FPRAS, in short) for approximating permanents of matrices with 
nonnegative entries^ [JSV^ . For an introduction to FPRAS, the reader may consult ^_H, p. xvii], [B and 
[JSVj . Unfortunately, the result of jJSV] is the best possible. If one could approximate by a FPRAS the 
permanent of matrices with (say) integer entries, then P = #P. 

There are however several approximation algorithms for general permanents, which do not yield a polynomial- 
time complexity bound. Perhaps the simplest one is based on the following theorem, noticed by several 
researchers (Hammersley, Heilman, Lieb, Gutman, Godsil). 

Theorem 2. Let A he a matrix and let B be the random matrix obtained from A by taking the square root of 
minimal argument of each non-zero entry and then multiplying each non-zero entry by an element o/{l, —1} 
chosen independently uniformly at random. Then E((det(i?))^) — per(v4). 

This leads to a Monte-Carlo algorithm for estimating the permanent. The algorithm was described and 
studied first by Karmarkar, Karp, Lipton, Lovasz and Luby (see [KKLLL] ) . Clearly, for general matrices the 
Monte-Carlo algorithm described above may have to run an exponential time; however, with regard to the 
connection to the quantum computing, the experimental study of the algorithm applied to the particular 
matrices which come from knot and link diagrams is an attractive task and it is our work in progress jointly 
with Petr Plechac. 

2. Proofs 

2.1. State sums and gadgets. Let us begin by recalling the state sum definition of the Jones polynomial 
of an oriented link, [J2l ITu| . Choose a diagram of a link L. Then, the Jones polynomial is given by the 
state sum 

J{L){q) = q-^'^'^D^)^^,oto(s)-rot,(s)^j^sgn(v)^^y 

S V 

where 

• A state s is a coloring of the edges of Z?l by or 1 such that around each vertex (positive or negative) 
oi Dl the coloring looks like one of the following possibilities: 



(1) 



(2) 



where edges colored by or 1 are depicted as dashed or solid respectively. 

The local weight R^{s) of a state s at a vertex v of (positive or negative) is given by 



+ 
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where q — q^^ ■ (These weights come from the P-matrix of the quantum group Uq{5i2) and we 
follow the conventions of [CP[ p.235]). 
• Lo{L) is the writhe oi Dl, that is the sum of the signs of the crossings oi D^. 
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• To define the rotation numbers roto(s) and roti(s) of a state s, we first define the the rotation number 
rot(a) of a curve a immersed in the plane to be the ip/{2'K), where ip is the total rotation angle of 
the tangent vector of a. (The direction of the counter-clockwise rotation is taken to be positive.) 
If s is a state of D^, then roti(s), for i — 1,2, is the sum of the rotation numbers of each of the 
z-colored circles in s. The rotation numbers can also be defined combinatorially as follows: smooth 
any 4-valent vertex of s as follows; 

The result is a collection of oriented planar circles, colored by or 1. roti(s) is the number of 
counter-clockwise 0-colored circles minus the number of clockwise i-colored circles. 

It will be convenient to include the phase factors q""^ in the _R-niatrix. This can be achieved as follows: 
Let Dl he the diagram obtained by isotoping Dl so that at each crossing the over and undercrossing arcs 
meet each other at a tangent where they cross. This means that the rotation number of a cycle c oi Dl 
equals to X^eec ^'*^^(^)' where rot(e) is the rotation number of the edge e of Dl corresponding to edge e of 
Dl. 

Now, for a state s and a vertex v of Dl, consider its two outgoing edges ei,e2 (see Figured]) and the 
rotation numbers ri = rot(ei) and r2 — rot(e2). Of course, ri and r2 depend on v. We now define the 
modified weights B^{s) of a state s at a vertex v by: 




(3) 



+ 
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With these weights, we have 

Lemma 2.1. Suppose that Dl is a diagram for a link L. Let Dl he the diagram obtained by isotoping Dl 
so that at each crossing the over and undercrossing arcs meet each other at a tangent where they cross. Then 

S V 

Proof. If s is a state of Dl and Si is the set of z-colored circles in s, i = 0, 1 then 

rot(e)(5e^s^ —Tot{e)5e,si 



^rotos-roti(s) _ ^rot(so)-rot(si) 



where the second equality follows since the four arcs entering an exiting a gadget are tangenial. This sum 
can be written as 

gEe starts at ^ rot(e)5e,so -mt{e)S^,s^ 



where (5e,c = 1 (respectively 0) if e lies in c (respectively does not lie in c). Thus, 

^roto(s)-roti(s) ^sgn(ti) ^^^^ ^ ^ J^^Eestarts at«™t(e)*e,so-™t(e)(5^,si^sgn(t))^^^ 



□ 

Note that in this proof we used the fact that in Dl , the over and undercrossing arcs meet each other at a 
tangent where they cross. We could avoid this tangent condition on the edges entering an exiting a crossing 
by adding factors to the the expressions B^{s) that are determined by the angles formed by the crossings 
in _Dl. See |Ba) for details on this type of construction. 
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Definition 2.2. Let Dl be a link diagram. Then we let Dl denote the directed, edge- weighted graph drawn 
in the plane which is obtained from Dl by replacing a neighbourhood of each crossing of Dl with the gadgets 
and their incident edges as shown in Figure [T] The edge weights W of Dl are given by the convention that 
all of the undecorated edges have weight 1, 

(a, b, c, d, e) = {{q - g-l), {q + q-^)/2, -q/2, q, 1/2), 

{v, w, X, y, z) = (1/2, - <?, (g + q-^)/2, -q-^/2, g-^), 
and each edge e coming from an edge e of Dl is assigned the weight g~2rot(e)^ 

Lemma 2.3. Let Dl be a link diagram of L and Dl be the graph constructed as described in Definition \2. 2i 
Then 

where V denotes the set of vertex- disjoint directed cycles that cover all vertices of Dl- 

Proof. Let s be a state of Dl and s be the corresponding state m Dl (the diagram constructed in Lemma [2.ip . 
There is a correspondence between the state 5 at a crossing v and sets of vertex-disjoint directed cycles that 
cover all vertices oi Dl- This correspondence is given in the following way: if a is an edge incident with v 
in Dl, and a is the corresponding edge in Dl, then the edge a is in a cycle if and only if a is colored by 1. 
The correspondence is given for a positive crossing in Appendix |^ From the appendix it is also readily seen 
that at a positive crossing 

B+{s)^ ^™t(ei)+rot(e,)) -Q^^^ 

p£'P{s,v) eSp 

where if V{s) denotes the set of vertex-disjoint directed cycles that cover all vertices of Dl that corresponds 
with the state s, then V{s,v) denotes the set of edges of V{s) that belong to the gadget at the crossing v 
of Dl- It is readily checked that the corresponding identity also holds at negative crossings. It then follows 
that 

p^V e^p s V 

which, by Lemma |2. 11 is equal to the Jones polynomial, as required. □ 

We make the following definition for the matrix M l that is used in Theorem [TJ 

Definition 2.4. Let Dl be a link diagram and Dl be the associated graph constructed as in Definition 1 2. 2 1 
Then we let Ml denote the adjacency matrix of Dl, whose («, j)-entry is the weight of the directed edge 

We can now prove Theorem [TJ 

Proof of Theorem]^ It is well known that if A is the 0, 1 adjacency matrix of a graph, then per(A) is the 
number of collections of vertex-disjoint directed cycles that cover all vertices of the graph. Similarly for 
general matrix A, per{A) = a(C), where the sum is over all collections of vertex-disjoint directed cycles 
that cover all vertices of the graph and a(C) — n{ij}GC^»i- Using the notation of Lemma [131 it follows 
that per(Mi) — J^pev Yieep The result then follows from Lemma[l?ni D 



Appendix A. Figures 



Case 1. q — ae + b 



Case 2. 1 1 




Case 3. 1 — bd + ac 




Case 4. q = d 




Case 5. q — q ^ = a 




+ 



Case 6. = de + c 
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